arXiv: 1509.05728vl [gr-qc] 18 Sep 2015 


KCL-PH-TH/2015-40 


Semiclassical solutions of generalized Wheeler-DeWitt cosmology 

Marco de Cesare^*, Maria Vittoria Gargiulo^^, Mairi Sakellariadou^^ 
^Department of Physics, King’s College London, 

University of London, Strand, London WC2R 2LS, U.K. and 
^Dipartimento di Fisica, I.N.F.N., Universitd di Salerno, I-SflOO Salerno, Italy 

We consider an extension of WDW minisuperpace cosmology with additional interaction terms 
that preserve the linear structure of the theory. General perturbative methods are developed and 
applied to known semiclassical solutions for a closed Universe filled with a massless scalar. The 
exact Feynman propagator of the free theory is derived by means of a conformal transformation 
in minisuperspace. As an example, a stochastic interaction term is considered and first order per¬ 
turbative corrections are computed. It is argued that such an interaction can be used to describe 
the interaction of the cosmological background with the microscopic d.o.f. of the gravitational field. 
A Helmoltz-like equation is considered for the case of interactions that do not depend on the in¬ 
ternal time and the corresponding Green’s kernel is obtained exactly.The possibility of linking this 
approach to fundamental theories of Quantum Gravity is investigated. 
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I. INTRODUCTION 

The Wheeler-DeWitt (WDW) equation represents the starting point of the canonical quantization program, also 
known as geometrodynamics. It was originally derived by applying Dirac’s quantization scheme to Einstein’s theory 
of gravity, formulated in Arnowitt-Deser-Misner (ADM) variables. It has the peculiar form of a timeless Schrodinger 
equation, whose solutions are functions defined on the space of three-dimensional geometries 


HtP = 0 . ( 1 ) 

The equation is a constraint imposed on physical states, encoding the invariance of the theory under reparametrization 
of the time variable. Together with the (spatial) diffeomorphism constraint, it expresses the principle of general 
covariance at the quantum level. This approach has however several limitations that prevented it from being accepted 
as a fully satisfactory theory of Quantum Gravity in its original formulation^. Nonetheless, some of the issues it 
raises, as for instance the problem of time {i.e. the lack of a universal parameter used to describe the evolution of the 
gravitational field), are actually problems encountered by all fundamental (nonperturbative) approaches to quantum 
gravity. Others stem instead from the dubious mathematical structure of the theory. Among the most essential ones 
in the latter category we mention the factor ordering problem in the Hamiltonian constraint and the construction of 
a physical Hilbert space. 

Despite its limitations, WDW was proved to be a valuable instrument in cases where the number of degrees of 
freedom is restricted a priori, i.e. in minisuperspace and midisuperspace models. In particular, its application to 
cosmological settings has been extremely useful to gain insight into some of the deep questions raised by a quantum 
theory of the gravitational held, such as the problem of time [2] and the occurrence of stable macroscopic branches 
for the state of the Universe (as in the consistent histories approach 0 ). 

One of the fundamental aspects of WDW is its linearity. This is indeed a consequence of the Dirac quantization 
procedure and a property of any hrst-quantized theory. One can nonetheless hnd in the literature several proposals for 
non-linear extensions of the geometrodynamics equation. They are in general motivated by an interpretation of if in 
Eq. § as a held operator (rather than a function) dehned on the space of geometries. This idea has been applied to 
Gosmology and found concrete realization in the old baby Universes approach [3]. More recently, a model motivated 
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by Group Field Theory (GFT) (see Ref. [S] for a recent outline) has been proposed in Ref. [5] that retains the same 
spirit, taking Loop Quantum Gosmology (LQG) [7] as a starting point^. 

It is worthwhile stressing that the WDW theory can be recovered in the continuum limit (see Ref. [S] for details) 
from a more fundamental theory such as LQG. In this sense it can be interpreted as an effective theory of Quantum 
Gosmology, valid at scales such that the fundamentally discrete structure of spacetime cannot be probed. Therefore 
WDW, far from being of mere historical relevance, is rather to be considered as an important tool to extract predictions 
from cosmological models when a continuum, semiclassical behaviour is to be expected. 

In what follows we consider as a starting point the Hamiltonian constraint of LQG, leading to the free evolution of 
the Universe, as discussed in the GFT-inspired model [^. Non-linear terms are allowed, which can be interpreted as 
interactions between disconnected homogeneous components of an isotropic Universe (scattering processes describe 
topology change). Instead of resorting to a third quantized formalism built on minisuperspace models, we follow a 
rather conservative approach which does not postulate the existence of Universes disconnected from ours. We therefore 
only take into account linear modifications of the theory, so as to guarantee that the superposition principle remains 
satisfied. Hence these extra terms can be interpreted as self-interactions of the Universe, as well as violations of the 
Hamiltonian constraint. 

We aim at studying the dynamics in a regime such that the free LQG dynamics (given by finite difference equation) 
can be approximately described in terms of differential equations. In this regime, the dynamics is given by a modified 
WDW equation. The additional interactions should be such that deviations from the Friedmann equation are small, 
and will therefore be treated as small perturbations. 

We consider a closed Friedmann-Lemaitre-Robertson-Walker (FLRW) Universe, for which semiclassical solutions are 
known explicitly |10j . Corrections to the solutions arising from the extra linear interactions are obtained perturbatively. 
Of particular interest is the case in which the perturbation is represented by white noise. In fact, this could be a way 
to model the effect of the discrete structure underlying spacetime on the evolution of the macroscopically relevant 
degrees of freedom. It should therefore be possible to make contact, at least qualitatively, between the usual quantum 
cosmology and the GFT cosmology, according to which the dynamics of the Universe is that of a condensate of 
elementary spacetime constituents. Considering a white noise term amounts to treating the vacuum fluctuations of 
the gravitational field as a stochastic process, an approach motivated by an analogy with Stochastic Electrodynamics 
(SED) [lIHIll- This will be done is Section [Vl} The analogy with SED has of course its limitations, since in the case 
at hand there is no knowledge about the dynamics of the vacuum, which should instead come from a full theory 
of Quantum Gravity. In spite of this limitation, the methods developed here are fully general, so as to allow for a 
perturbative analysis of the solutions of the linearly modified WDW equation for any possible form of the additional 
interactions. 

In the present work we restrict to a closed Universe, for which wave packet solutions were constructed in Ref. [10] . 
In the same work, it was shown that it is possible to construct a quantum state whose evolution mimics a solution 
of the classical Friedmann equation, and is given by a quantum superposition of two Gaussian wave packets (one for 
each of the two phases of the Universe, expanding and contracting) centered on the classical trajectory. We build on 
the results of Ref. [10) . considering the new interaction term as a perturbation, and thus determining the corrections 
to the motion of the wave packets. 

The rest of the paper is organized as follows: In Sec. H, we show how our model is motivated by LQG and GFT 
inspired cosmology [6] in the continuum limit (large volumes). In Sec. HI, we review the construction of the wave 
packet solutions done in Ref. m- These solutions will be used in the subsequent sections as the unperturbed states 
describing the evolution of a semiclassical Universe. In Sec. IV, using the scalar field to define an internal time, we 
develop a framework for time-independent perturbation theory. The inverse of the Helmoltz operator corresponding to 
the free theory is computed exactly and turns out to depend on a real parameter, linked to the choice of the boundary 
conditions at the singularity. An analogue of Ehrenfest’s theorem for the evolution equation of the expectation values 
of observables is then given. In Sec. V, the Eeynman propagator of the WDW free field operator is evaluated exactly. 
This is accomplished by using a conformal map in minisuperspace, which reduces the problem to that of finding the 
Klein-Gordon propagator in a planar region with a boundary. In Sec. VI, we consider the case in which the additional 
interaction is given by white noise. In Sec. VH, we discuss the role of the Noether charge and the choice of the inner 
product. Finally, we review our results and their physical implications in Sec. VIII. 


^ Yet another interpretation of the dynamical equation on minisuperspace has been recently advocated in Ref. According to this 
point of view, it should rather be interpreted as some analogue of a (non-linear) hydrodynamical equation, for which the superposition 
principle does not hold in general. 
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II. FROM THE LQC FREE THEORY TO WDW 

Let us briefly show how the WDW equation is recovered from the Hamiltonian constraint of LQC following Ref. 

For our purposes it is convenient to consider a massless scalar field (j), minimally coupled to the gravitational field. 
This choice has the advantage of allowing for a straightforward deparametrization of the theory, thus defining a clock. 
The Hamiltonian constraint has the general structure ElISl 

JC'tlj{i^,(t)) = -B{iy){Q + dl)tp{v,(l))=0, ( 2 ) 

where ^ is a wave function on configuration space, 0 is a finite difference operator acting on the gravitational sector 
in the kinetic Hilbert space of the theory and v denotes the generic eigenvalue of the volume operator. The 

discreteness introduced by the LQC formulation does not affect the matter sector, which is still the same as in 
the continuum WDW quantum theory. The gravitational sector of the Hamiltonian constraint operator of LQC in 
“improved dynamics” reads 

- B{v)Qij){v,(j)) = A{v)'il){v + vo,(t>) + C{v)il}{v,(j)) + D{v)ij){v - i^o,4>) , (3) 

where the finite increment z/q represents an elementary volume unit and A^B^C,D are functions which depend on 
the chosen quantization scheme. In order to guarantee that 0 is symmetric® in z/, the coefficients must satisfy the 
D{v) = A(y — i/q) condition [^; it holds in both the k = 0 and the k = 1 case. 

It is a general result of LQC that WDW can be recovered in the continuum (be. large volume) limit [TS]. In 
particular, it was shown in Ref. [5] that for k = 1 one recovers the Hamiltonian constraint of Ref. m- In fact it turns 
out that 0 can be expressed as the sum of the operator (0o) relative to the k = 0 case and a independent potential 
term (be. diagonal in the v basis) as 


0 — 00 + 


,4/3 

3^4/3 I I 


(4) 


In the above expression K = 


2\/2 


, 7 is the Barbero-Immirzi parameter of Loop Quantum Gravity (LQG), G the 


gravitational constant and Iq is the cube root of the fiducial volume of the fiducial cell on the spatial manifold in the 
k = 1 model. The latter can be formally sent to zero in order to recover the k = 0 case. 

Restricting to wave functions which are smooth and slowly varying in z/, we obtain the WDW limit of the 
Hamiltonian constraint 


0o'0(g<(>) = -l2nG{L'd^f'il;{iy,cj)) , (5) 

which is exactly the same constraint one obtains in WDW theory. Thus, LQC naturally recovers the factor ordering 
(also called factor ordering, in the sense that the quantum constraint operator is of the form where 

G"^^ is the inverse WDW metric and denotes the covariant derivative associated with Gab) which was obtained 
in Ref. [2] under the requirement of field reparametrization invariance of the minisuperspace path-integral. Since z/ 
represents a proper volume, it is proportional to the volume of a comoving cell with linear dimension equal to the 
scale factor a 


V oc. . 


Introducing the variable a = logo, we rewrite the constraint operator K, as 


t = e-®“ 




( 6 ) 

(7) 


® The actual way in which WDW represents a large volume limit of LQC is put in clear mathematical terms in Refs. 11511161 . where the 
analysis is based on a special class of solvable models (sLQC). 

^ In the framework of LQC, the Hamiltonian constraint contains the gravitational connection c. However, only the holonomies of the 
connections are well defined operators, hence to quantize the theory we replace c by sin/ic//2, where Ji represents the “length” of the 
line segment along which the holonomy is evaluated. Originally Ji was set to a constant /iq? related to the area-gap. To cure severe issues 
in the ultraviolet and infrared regimes which plague the /io quantization, a new scheme called “improved dynamics” was proposed El- 
In the latter, the dimensionless length of the smallest plaquette is fi. 

^ In the large volume limit it will be formally self-adjoint w.r.t. the measure B{u)du. 
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Let us consider a modified dynamical equation of the form 


(/)) 


E 


d(j)' giy, (j), 4>') = 0 


which, in the continuum limit v leads to 


da^ 


d(jp 



da' d(p' 


diy 

da 



We will assume that the interaction g is local in minisuperspace, i.e. g{v, v'; cj), (j)') = g(i/, (fjSiv 
the properties of the Dirac delta function 


( 8 ) 

(9) 

vi)5{(j) — (j)'). Given 


5{v — v') 



( 10 ) 


Equation ([^ reduces to a Klein-Gordon equation with space and time dependent potential e^°‘g{a, (j)); note that g is 
kept completely general. 

Since the interaction term represented by the r.h.s. of Eq.(|^ is unknown, we cannot determine an exact solution of 
the equation without resorting to a case by case analysis. However, since the solutions of the WDW equation in the 
absence of a potential are known explicitly, we will adopt a perturbative approach. The method we develop is fully 
general and can thus be applied for any possible choice of the function g. 

We formally expand the wave function and the WDW operator in terms of a dimensionless parameter A (that serves 
book-keeping purposes and will be eventually set equal to 1). 


with the definitions 


We therefore have 


jP = jpiO) y + . . . , 
and f = to + Ae3“ J g , 

f = i — - — - 

°“^aa2 502 ® y- 


= 0 , 


( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 


The zero-th order term is a solution of the wave equation with an exponential potential, which was obtained 
in Ref. |10j and will be reviewed in the next section. If we were able to invert Tg, we would get the wave function 
corrected to first order. However, finding the Green’s function is not straightforward in this case as it would be for 
k — 0 (where the kinetic operator is just the d’Alembertian, whose Green’s kernels are well known for all possible 
choices of boundary conditions). Moreover, as for the d’Alembertian, the Green’s kernel will depend on the boundary 
conditions. The problem of determining which set of boundary conditions is more appropriate depends on the physical 
situation we have in mind and will be dealt with in the next sections. 

At this point, we would like to point out the relation between our approach and the model in Ref. [6], where a third 
quantization perspective is assumed. The action on minisuperspace that was considered in Ref. [B] reads® 

S'[V’] = 5'free[0] ^ ^ /d^i . . .d0„ /("^(^'*,0i)]j0>(z/j,0j) , (16) 

Tl. J 

n U\..Mn J 


where the first term gives the dynamics of the free theory, namely a homogeneous and isotropic gravitational back¬ 
ground coupled to a massless scalar field 


s„„w.] = y: / d0 '0(j^, 0) ; 


(17) 


® The model was originally formulated for k = 0, but it admits a straightforward generalization to include the case k = 1. 
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1C is the Hamiltonian constraint in LQC and the terms containing the functions represent additional interactions 
that violate the constraint. This is a toy model for Group Field Cosmology [6], given by a GFT with group G = U{1) 
and v a Lie algebra element. The free dynamics depends on the specific LQC model adopted. However, the continuum 
limit should be the same regardless of the model considered and must give the WDW equation for the corresponding 
three-space topology. The WDW approach to quantum cosmology should therefore be interpreted as an effective 
theory, valid at scales such that the discreteness introduced by the polymer quantization cannot be probed. We will 
show that, even from this more limited perspective, the action Eq. (16) leads to novel effective theory of quantum 
cosmology that represent modifications of WDW. 

It was hinted in Ref. that the additional interactions could also be interpreted as interactions occurring between 
homogeneous patches of an inhomogeneous Universe. Another possible interpretation is that they actually represent 
interactions among different, separate, Universes. The latter turns out to be a natural option in the framework of 
third quantization (see Ref. m and references therein), which naturally allows for topology change^. 

We restrict our attention to the quadratic term, that can be interpreted as a self interaction of the Universe. This 
is in fact worth considering even without resorting to a third-quantized cosmology and can in principle be generalized 
to include non-locality in minisuperspace (this situation will not be dealt with in the present work). 


III. ANALYSIS OF THE UNPERTURBED CASE 


We review here the construction of wave packet solutions for the cosmological background presented in Ref. m- 
The Wheeler-DeWitt equation for a homogeneous and isotropic Universe (compact spatial topology, fc = 1) with a 
massless scalar held is 


/ 

\ da^ dcjP' 


Aa 


i’AA) = 0 ■ 


We impose the boundary condition 


lim (/)) = 0 , 

Oi—^OCi 


(18) 


(19) 


necessary in order to reconstruct semiclassical states describing the dynamics of a closed Universe, since regions of 
minisuperspace corresponding to arbitrary large scale factors are not accessible. 

Equation (191 can be solved by separation of variables 


4>) = Nk Ck{a)ipk{(t>) , 


leading to 


52 

+ k (fik = 0 , 


52 


^Cfc-(e4“-fc2)c,=0. 


( 20 ) 

( 21 ) 

( 22 ) 


Note that ipk{,o.A) stands for the solution of Eq. (18) and should not be confused with the Eourier transform of ip, 
for which we will use instead the notation ijj(a, k). In Eq. (20) above, Nk is a normalization factor that depends on k, 
whose value will be hxed later. 


Let us proceed with the solutions of Eqs. (21), (|22[). Equation (21) yields complex exponentials as solutions 


^k = 


(23) 


while Eq. (22) has the same form as the stationary Schrodinger equation for a non-relativistic particle in one dimension, 
with potential V{a) = — k^ and zero energy. One then easily remarks that the particle is free for a —> — 00 , 

whereas the potential barrier becomes inhnitely steep as a takes increasingly large positive values. Given the imposed 
boundary condition, Eq. (22) admits as an exact solution the modihed Bessel function of the second kind (also known 
as MacDonald function) 


Gk{A — kCik/2 


( 24 ) 


^ For an example of topology change in the old “Baby Universes” literature see e.g. Ref. [4]. 
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Wave packets are then constructed as linear superpositions of the (appropriately normalized) solutions 


tjj{a,(j))= / Ak'ipk{a,(j))A{k) , 


(25) 


with an appropriately chosen amplitude A(k ). Since the Ck{a) are improper eigenfunctions of the (one-parameter 
family of) Hamitonian operator(s) in Eq. (22), they do not belong to the space of square integrable functions on the 
real line L^(IR). However, it is still possible to define some sort of normalization by fixing the oscillation amplitude of 
the improper eigenfunctions for a —>■ —oo. For this purpose, we recall the WKB expansion of Eq. (24) 




o4a'1-1/4 


COS 


f ^ arccosh -^ ~ j 

V2 e2“ 2 4/ ’ 


a very accurate approximation for large values of k. We therefore set 




(26) 


(27) 


which for large enough values of k gives elementary waves with the same amplitude to the left of the potential barrier. 
Note that for small k, the amplitude A{k) will still exhibit dependence on k. 

Let us assume a Gaussian profile for the amplitude A{k) in Eq. (25) 


A{k) = 


1 


(fc-fc)^ 




(28) 


where k should be taken large enough so as to guarantee the normalization of the function ipkio:, 4>). One then finds that 
the solution has the profile shown in the Fig.[2 The solution represents a wave packet that starts propagating from a 
region where the potential vanishes {i.e. at the initial singularity) towards the potential barrier located approximately 
at = I log fc, from where it is reflected back. For such values of k the wave packet is practically completely reflected 
back from the barrier. The parameter b gives a measure of the semiclassicality of the state, i.e. it accounts for how 
much it peaks on the classical trajectory. We remark that the peak of the wave packet follows closely the classical 
trajectory, 


and refer the reader to Fig. 


k 

cosh(2(/)) 


(29) 


IV. TIME INDEPENDENT PERTURBATION POTENTIAL 


In the case in which the potential in the r.h.s of Eq. does not depend on the scalar field (that we interpret as an 
internal time) we can resort to time-independent perturbation theory to calc ulat e the corrections to the wave function. 
Namely, we consider the representation of the operator Tq defined in Eq. (13) in Fourier space, or equivalently its 


action on a monochromatic wave, so that we are lead to the Helmoltz equation in the presence of a potential 


(30) 


where the potential cannot be considered as a small perturbation with respect to the standard Helmoltz equation. 
In fact, reflection from an infinite potential barrier requires boundary conditions that are incompatible with those 
adopted in the free particle case. We therefore have to solve Eq. (301 exactly. 


In order to compute the first perturbative corrections, we need to solve Eq. (15), which in Fourier space leads to 

( k) = Uk) , (31) 


upon defining U{a) = —e^^g{a). Equation (31) is easily solved once the Green’s function of the Helmoltz operator 


on the l.h.s is known. The equation for the Green’s function is 


(— -k - e'‘“ 


^ G'^{a,a')=5{a-a') . 


(32) 
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Figure 1: Absolute square of the “wave function” of 
the Universe corresponding to the choice of 
parameters b = 1, k = 10. Lighter shades correspond 
to larger values of the wave function. During 
expansion and recollapse the evolution of the 
Universe can be seen as a freely propagating wave 
packet. From the plot it is also evident the reflection 
against the potential barrier at ^ log k, where 

the wave function exhibits a sharp peak. 



Figure 2: Classical trajectory of the Universe in 
minisuperspace. It is closely followed by the peak of 
semiclassical states, as we can see by comparison 
with Fig. 


The homogeneous equation admits two linearly independent solutions, which we can use to form two distinct linear 
combinations that satisfy the boundary conditions at the two extrema of the interval of the real axis we are considering. 
The remaining free parameters are then fixed by requiring continuity of the function at a' and the condition on the 
discontinuity of the first derivative at the same point. 

Equation (301 has two linearly independent solutions p{a) and h(a) given by 


and 


We thus make the following ansatz 


P{a) =K^k{ — 




(33) 


(34) 


G^{a,a') = 


lP{a) 


a> a' 


5 h(a) + 77 p[a) a < a' 


(35) 


and note that G^{a,a’) satisfies the boundary condition Eq. (191 by construction. 

Since we do not know what is the boundary condition for a —>■ — 00 , i.e. near the classical singularity®, we will 
not be able to fix the values of all constants 7,(5, ??. We will therefore end up with a one-parameter family of Green’s 
functions. 


In fact, even popular choices like the Hartle-Hawking no-boundary proposal |19II20| or the tunneling condition proposed by Vilenkin ED, 
only apply to massive scalar fields. 











The Green’s function must be continuous at the point a = a', hence 

7 p{a') = (5 h(a') + rj p{a') . (36) 

Moreover, in order for its second derivative to be a Dirac delta functional, the following condition on the discontinuity 
of the first derivative must be satisfied 


7 p'{a') — 6 h'(a') — rj p'{a') = 1 . 


(37) 


Upon introducing a new constant U = 7 — 77 , we can rewrite Eqs. (36), (37) in the form of a Kramer’s system 

f D p{a') — S h{a') = 0 , 

I U p'{a') — 6 h'{a') = 1 , 


which admits a unique solution, given by 


D = - 
5 = - 


h{a') 
W{a') ’ 

p{a') 


W{a') ’ 

where W{a') = p{a')h'{a') — h{a')p'{a') is the Wronskian. Hence the ansatz (35) is rewritten as 


(38) 

(39) 


G'=’’'(a,a') = 


p{a) {p 




W{a') 


a> a' 


where we have explicitly introduced the parameter p in our notation for the Green’s function in order to stress its 
non-uniqueness. Note that the phase shift at —00 varies with p. 

Hence the solution of Eq. (31) reads 




(40) 


The possibility of studying the perturbative corrections using a decomposition in monochromatic components is viable 
because of the validity of the superposition principle. This is in fact preserved by additional interactions of the type 
considered here, which violate the constraint while preserving the linearity of the wave equation. Note that in general, 
modifications of the scalar constraint in the WDW theory would be non-linear and non-local, as for instance within 
the proposal of Ref. [5], where classical geometrodynamics arises as the hydrodynamics limit of GET. However, such 
non-linearities would spoil the superposition principle, hence making the analysis of the solutions much more involved. 
The time dependence of the perturbative corrections is recovered by means of the inverse Fourier transform of 

Eq. (Idol) 


i/,(i)(a,(/,) = y ^^W(^a,k)e 


(41) 


Expectation values of observables can be dehned using the measure determined by the time component of the 
Noether current (whose dehnition will be given in Sec. VH) as 


(/) = y ■ 


Using the conservation of the Noether current we can then derive an analogue of Ehrenfest theorem, namely 

d 


dcj) 


if) = J da + fd^j^) = J da + fd^ja) ■ 


(42) 


(43) 


When the observable / does not depend explicitly on the internal time (j), the first term in the integrand vanishes. 
Gonsidering for instance the scale factor a = e“ then, after integrating by parts, we get 
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In an analogous fashion, one can show that 


d(/) 



da . 


(45) 


Formulae (43) and (45), besides their simplicity, turn out quite handy for numerical computations, especially when 
dealing with time independent observables. In fact they can be used to compute time derivatives of the averaged 
observables without the need for a high resolution on the (f) axis, i.e. they can be calculated using data on a single 
time-slice. Expectation values can therefore be propagated forwards or backwards in time by solving first order 
ordinary differential equations. 


V. FIRST PERTURBATIVE CORRECTION FOR TIME-DEPENDENT POTENTIALS 


In the previous section we considered a time-independent perturbation, which can be dealt with using the Helmoltz 
equation. This is in general not possible when the perturbation depends on the internal time. In order to study the 
more general case we have to resort to different techniques to find the exact Green’s function of the operator Tg. Since 
the potential 6"^“ breaks translational symmetry, Fourier analysis, which makes the determination of the propagator 
so straight-forward in the k = 0 case (where the Hamiltonian constraint leads to the wave equation), is of no help. 
Let us perform the following change of variables 

X = cosh(2(/)) , Y = sinh(2(/)) , (46) 

which represents a mapping of minisuperspace into the wedge X > |V|. The minisuperspace interval (corresponding 
to DeWitt’s supermetric) can be expressed in the new coordinates as 

d02_da2 = _^(dr2_d^2) ^ (47) 

with (X^ — Y^)~^ the conformal factor®. As an immediate consequence of conformal invariance, uniformly expanding 
(contracting) Universes are given by straight lines parallel to X — Y {X = —Y) within the wedge. The “horizons” 
X = — Y and X = Y represent the initial and final singularity, respectively. It is worth pointing out that the classical 
trajectory Eq. ( |29[ ) takes now the much simpler expression 

k = 2X , (48) 


i.e. classical trajectories are represented by straight lines parallel to the Y axis and with the extrema on the two 
singularities. 


In the new coordinates X, Y the operator Tq in Eq. (13) reads 

To = (A® - y2)(ai - - 4) , 


(49) 


which is, up to the inverse of the conformal factor, a Klein-Gordon operator with = 4. This is a first step towards 
a perturbative solution of Eq. (15), which we rewrite below for convenience of the reader in the form 


where 


= C/(a,(^)V'(o) , 

U{a,<p) = -e^°‘g{a,(j)) . 


(50) 

(51) 


Note that the above is the same equation as the one considered in the previous section, but we are now allowing for 
the interaction potential to depend also on (j). Given Eq. (49), we recast Eq. (50) in the form that will be used for the 
applications of the next section, namely 


(dx -dy- 


U{X,Y) 
(X2 - r2) 




( 52 ) 
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Recall that any two metrics on a two-dimensional manifold are related by a conformal transformation. 
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The formal solution to Eq. (52) is given by a convolution of the r.h.s with the Green’s function satisfying suitably 
chosen boundary conditions 


V'(i) (X, Y)=J dX'dF' GiX, Y; X', (X', Y') . (53) 

The Green’s function of the Klein-Gordon operator in free space is well-known for any dimension D (see e.g. Ref. |22)i. 
In = 2 it is formally given by^° 


G{X,Y-X',Y') = I 


(27r)2 (fcf. - fc|) - 4 


and satisfies the equation 


(5i -d^- 4)G(X, y; X', Y') = 5{X - X')5{Y - Y') . 


(54) 


(55) 


Evaluating (54) explicitly using Feynman’s integration contour, which is a preferred choice in the context of a third 
quantization, we get [HI |H] 


G{X,Y;X',Y') = -i0(s)iff (2^5) - ±g{-s)Ko{2y^s) , (56) 

where we introduced the notation s = (F — F')^ — {X — X'Y for the interval. However, the present situation 
is distinguished from the free case, since there is a physical boundary represented by the edges of the wedge. The 
boundary conditions must be therefore appropriately discussed. A preferred choice is the one that leads to the Feynman 
boundary conditions in the physical coordinates In the following, we will see the form that these conditions 

take in the new coordinate system, finding the transformation laws of the operators that annihilate progressive and 
regressive waves. 

We begin by noticing that the generator of dilations in the (A, F) plane acts as a tangential derivative along the 
edges 


X dx + F dy = dt . 


(57) 


Moreover on the upper edge X = Y (corresponding to the big crunch) we have 


dt\x^Y = ^ {da + d^) , 

while on the lower edge (corresponding to the initial singularity) we have 

dt\x=-Y = 4 


(58) 


Therefore the boundary conditions 


dt\x=YG = 9 , 


't\X=-Y 


d^) . 

(59) 

G = 0 

(60) 


are equivalent to the statement that the Green’s function is a positive (negative) frequency solution of the wave 
equation at the final (initial) singularity. This is in agreement with the Feynman prescription and with the fact that 
the potential e“ vanishes at the singularity a —> — oo. 

There is a striking analogy with the classical electrostatics problem of determining the potential generated by a 
point charge inside a wedge formed by two conducting plates (in fact it is well-known that the electric field is normal 
to the surface of a conductor, so that the tangential derivative of the potential vanishes). The similarity goes beyond 
the boundary conditions and holds also at the level of the dynamical equation. In fact, after performing a Wick 
rotation Y ^ —iY Eq. (55) becomes the Laplace equation with a constant mass term, while the operator dt defined 


in Eq. (57) keeps its form. 


After performing the Wick rotation the problem can be solved with the method of images. Given a source (charge) 
at point Pq = {r',9'), three image charges as in Fig. are needed to guarantee that the boundary conditions are met. 


Notice that here Y plays the role of time. 
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Figure 3: The whole of minisuperspace is 
conformally mapped into the wedge. The straight 
lines X = ±Y correspond to the two horizons, while 
the vertical line is the classical trajectory of a closed 
Universe. 



Figure 4: The Green’s function is a function of the 
point Q and the source Pq. The positions and the 
charges of the images Pj ’s are determined so as to 
satisfy the boundary conditions. 


The Euclidean Green’s function as a function of the point Q = (r, 9) and source Pq (hence m = 2) is shown to be 
given by 


G{r,9,r',9') =-^ (Ko{mPoQ) - Ko{mPiQ) + KoimP 2 Q) - Ko^mPsQ)) . (61) 

The quantities PjQ represent the Euclidean distances between the charges and the point Q. The Lorentzian Green’s 
function is then recovered by Wick rotating all the Y time coordinates, i.e. those of Q and of the Pj’s. In order to 
obtain this solution we treated the two edges symmetrically, thus maintaining time reversal symmetry. The Green’s 
function with a source Pq within the wedge vanishes when Q is on either of the two edges. In fact, this can be seen 
as a more satisfactory way of realizing DeWitt’s boundary condition, regarding it as a property of correlators rather 
than of states^^. 


VI. TREATING THE INTERACTION AS WHITE NOISE 


The methods developed in the previous sections are completely general and can be applied to any choice of the 
extra interaction terms using Eq. (53). In this section we will consider, as a particularly simple example, the case in 


which the additional interaction is given by white noise. Besides the mathematical simplicity there are also physical 
motivations for doing so. In fact, we can make the assumption that the additional terms that violate the Hamiltonian 
constraint can be used to model the effect of the underlying discreteness of spacetime on the evolution of the Universe. 

Some approaches to Quantum Gravity (for instance GFT) suggest that an appropriate description of the gravita¬ 
tional interaction at a fundamental level is to be given in the language of third quantization |26j . At this stage, one 
might make an analogy with the derivation of the Lamb shift in Quantum Electro-Dynamics (QED) and in effective 
stochastic approaches. As it is well known, the effect stems from the second quantized nature of the electromagnetic 
field. However, the same prediction can also be obtained if one holds the electromagnetic field as classical, but impose 
ad hoc conditions on the statistical distribution of its modes, which necessarily has to be the same as that corre¬ 
sponding to the vacuum state of QED. In this way, the instability of excited energy levels in atoms and the Lamb 


DeWitt originally proposed the vanishing of the wave function of the Universe at singular metrics on superspace, suggesting in this 
way that the singularity problem would be solved a priori with an appropriate choice of the boundary conditions. However, there are 
cases (see discussion in Ref. mus]) where DeWitt’s proposal does not lead to a well posed boundary value problem and actually 
overconstrains the dynamics. For instance, the solution given in Ref. m that we discussed in Section [lll| satisfies it for a —>■ oo but not 
at the initial singularity, since all the elementary solutions C'fc(Q) are indefinitely oscillating in that region. Imposing the same condition 
on the Green’s function does not seem to lead to such difficulties. 
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shift are seen as a result of the interaction of the electron with a stochastic background electromagnetic held [HIS]. 
However, in the present context we will not resort to a third quantization of the gravitational held, but instead put 
in ad hoc stochastic terms arising from interactions with degrees of freedom other than the scale factor. Certainly, 
our position here is weaker than that of Stochastic Electro-Dynamics (SED) (see Ref. [T^l E] and the works cited 
above), since a fundamental third quantized theory of gravity has not yet been developed to such an extent so as to 
make observable predictions in Cosmology. Therefore we are unable to give details about the statistical distribution 
of the gravitational degrees of freedom in what would correspond to the vacuum state. Our model should henceforth 
be considered as purely phenomenological and its link to the full theory will be clarihed only when the construction 
of the latter will eventually be completed. 

To be more specific, we treat the function g in the perturbation as white noise. Stochastic noise is used to describe 
the interaction of a system with other degrees of freedom regarded as an environment^'^. Hence we write 

{9{a,(t>)) = 0 , 

( 5 ( 0 :, <(>')) = <^(a - a')5{(j) - cf)') , 

where (_) denotes an ensemble average. It is straightforward to see that 

= 0 , (64) 


(62) 

(63) 


which means that the ensemble average of the corrections to the wave function vanishes. 

A more interesting quantity is represented by the second moment 

J^{Xi,Yi,X2,Y2) = ((V’(i)(Xi,yi))*V'^i)(X2,F2)) . (65) 


In fact when evaluated at the same two points J-{X,Y) = J-{X, Y ; X, Y), it represents the variance of the statistical 
fluctuations of the wave function at the point {X,Y). Using Eqs. (15), (63) and (65) we get 


J^(X,r) = 16 [ dX'dY' {X'"^-Y'‘^)^\jP<^°\X',Y')\‘^\G{X - X',Y-Y')\‘^ . 
Jx>\Y\ 


( 66 ) 


In order to obtain the correct expression of the integrand, one needs to use the transformation properties of the Dirac 
distribution, yielding 


{g{a, 4>)g{a', 4>')) = 2(^2 - Y^)5{X - X')5{Y - Y') . 


(67) 


Notice the resemblance of Eq. (65) with the two-point function evaluated to first order using Feynman rules for a 
scalar held in two dimensions interacting with a potential (X^ — Y ^)2 Following this analogy, the variance X 

can be seen as a vacuum bubble. 

Equation (66) implies that a white noise interaction is such that the different contributions to the modulus square 
of the perturbations add up incoherently. 


VII. CONSERVED CHARGE AND THE CHOICE OF THE INNER PRODUCT 


Let us briefly discuss the role of the Noether charge and the choice of the inner product. Equation (18) is just a 
Klein-Gordon equation (with ip interpreted as time^^) with a potential V{a) = e'^°‘ — 

Equation (18) can be seen as the Euler-Lagrange equation of the Lagrangian 


dip* dip 
da da 


dip* dip 
d(p dcp 


+ ■ 


( 68 ) 


For a derivation of the Schrodinger-Langevin equation using the methods of stochastic quantization we refer the reader to Ref. m 
We are not aware of other existing works suggesting the application of stochastic methods to Quantum Cosmology. There are however 
applications to the fields of Classical Cosmology and inflation and the interested reader is referred to the published literature. 

We point out that another reason to prefer (/> as a time variable instead of a is that this choice leads to a potential that is bounded from 
below. Even from a classical perspective a is not a good choice as a global time since it is not a monotonic function of the proper time. 
However, it can still be used as a local time (a consistent way of using local times was studied in Ref. ([28]). In fact this is a problem 
one has to deal with for general, non-deparameterizable systems. 
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Since the Lagrangian is manifestly C/(l) invariant, the following conservation law holds as a consequence of Noether’s 
theorem 


= 0, 


where^^ 


i/x = -i 




The corresponding Noether charge (obtained following the standard procedure) reads 



(69) 


(70) 


(71) 


and is conserved under time evolution 


Clearly, the conservation of Q is in general incompatible with that of the norm of ^/>, defined by 

poo 

II^IIl2(r) = / da IV'P . 


(72) 


(73) 


Furthermore, the Noether charge Q has no definite sign unless one considers only solutions with either positive or 
negative frequency, since a generic state admits a decomposition in terms of positive and negative frequency solutions 
of the free WDW equation. Such solutions span superselection sectors that are preserved by the Dirac observables 
identified in Ref. [55], and states in a given superselection sector satisfy a Schrodinger equation with a Hamiltonian 
i-v/©- Within each sector, the conservation of the Hilbert norm holds true and Eq. (721 is an obvious consequence 


of energy conservation. As far as the free dynamics is concerned, it is perfectly legitimate to work in a superselection 
sector and use it to construct the (kinematical) Hilbert space^® |9l [29]. However, when the theory includes either 
non-linear interactions or potentials that depend on (internal) time, superselection no longer holds and both sectors 
have to be taken into account, thus naturally leading to a third quantization |18j . The non-equivalence of Q and 
can be easily seen for a wave packet (25), in which case one has 


Q = 


dk 2k\A{k)\‘^ 


(74) 


Note that when A{k) is supported only on one side of the real line, Eq. (72) can be interpreted as a conservation 
law for the average momentum^®. A similar interpretation is not available instead when both positive frequency and 
negative frequency solutions are present in the decomposition of the wave packet. Eor the reasons given above, a naive 
probabilistic interpretation of the theory along the lines of ordinary Quantum Mechanics does not seem possible in 
general. It is thus preferable to avoid using the designation of wave function for the classical field ip. 


The Noether charge of the semiclassical Universe considered in Section HI can actually be computed analytically, 
and since it is conserved one can evaluate it in the simplest case, i.e. when the potential vanishes for a —>■ —oo. In 
fact when both a and (j) are in a neighborhood of —oo, i.e. close to the initial singularity, one can approximate the 


We introduced a compact notation to denote the antisymmetric combination of right and left derivatives, namely 

9 ^ 8 _ 8 
8(f> 8(t> d4> 

The relation between different choices for the physical inner product and the Green’s functions in LQC was explored in Ref. m- 
In the QFT context the same relation is interpreted as a relativistic normalization, which amounts to having a density of 2k particles 
per unit volume. 
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exact solution with a Gaussian wave packet (given by the WKB approximation, see Appendix]^, which describes 
the evolution of the Universe in the expanding phase. We obtain 


V' ~ cr&U - exp < - 


k[ 2 (j) + arccosh { ^ ] + \l ^ ^ - 7 ^ 


exp 


' 2 


^ arccosh ( 

^ 20 : 


(75) 


and then using Eq. (70) the charge density reads 


U = |c^| 7r6"^fcexp 


( 20 + arccosh ( ^ 


(76) 


For large negative values of a one can approximate the inverse hyperbolic function in the argument of the exponential 
as 


arccosh ~ log(2fc) — 2a , (77) 

and hence reduce the evaluation of the charge to that of a Gaussian integral, leading to the result 

Q = \c^h'K^I'^k . (78) 

From the above analysis we conclude that the quantum conservation law leads, in the classical limit, to momentum 
conservation. We recall that classically the momentum is related to the derivative of the scalar field w.r.t. proper time, 
as 


k = a^<p . (79) 

VIII. CONCLUSIONS 

Motivated by LQC and GFT we considered an extension of WDW minisuperspace cosmology with additional 
interaction terms representing a self-interaction of the Universe. Such terms can be seen as a particular case of the 
model considered in Ref. 0 , which was proposed as an approach to quantum dynamics of inhomogeneous cosmology. 
In general, inhomogeneities would lead to non-linear differential equations for the quantum field. This is indeed the 
case when the “wave function” of the Universe is interpreted as a quantum field or even as a classical field describing 
the hydrodynamics limit of GFT. 

In the framework of a first quantized cosmology, we only considered linear modifications of the theory in order 
to secure the validity of the superposition principle. Our work represents a first step towards a more general study, 
that should take into account non-linear and possibly non-local interactions in minisuperspace. However, the way such 
terms arise, their exact form, and even the precise way in which minisuperspace dynamics is derived from fundamental 
theories of Quantum Gravity should be dictated from the full theory itself. 

Assuming that the additional interactions are such that deviations from the Friedmann equation are small, we 
developed general perturbative methods which allowed us to solve the modified WDW equation. We considered a 
closed FLRW Universe filled with a massive scalar field to define an internal time, for which wave packets solutions 
are known explicitly and propagate with no dispersion. A modified WDW equation is then obtained in the large volume 
limit of a particular GFT inspired extension of LQC for a closed FLRW Universe. Perturbative methods are then used 
to find the corrections given by self-interactions of the Universe to the exact solution given in Ref. m- To this end, the 
Feynman propagator of the WDW equation is evaluated exactly by means of a conformal map in minisuperspace and 
(after a Wick rotation) using the method of the image charges that is familiar from electrostatics. This is potentially 
interesting as a basic building block for any future perturbative analysis of non-linear minisuperspace dynamics for a 
closed Universe. 

A Helmoltz-like equation was obtained from WDW when the extra interaction does not depend on the internal 
time. Its Green’s kernel was evaluated exactly and turned out to depend on a free parameter 77 related to the choice 
of boundary conditions. Further research must be carried over to link 77 to different boundary proposals. 

We illustrated our perturbative approach in the simple and physically motivated case in which the perturbation is 
presented by white noise. In this phenomenological model the stochastic interaction term can be seen as describing 
the interaction of the cosmological background with other degrees of freedom of the gravitational field. Calculating 
the variance of the statistical fluctuations of the wave function, we found that a white noise interaction is such that 
the different contributions to the modulus square of the perturbations add up incoherently. 
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A. WKB APPROXIMATION OF THE ELEMENTARY SOLUTIONS AND ASYMPTOTICS OF THE 

WAVE PACKETS 


Since Eq. (22) has the form of a time-independent Schfodinger equation, it is possible to construct approximate 


solutions using the WKB method. 

For a given A:, we divide the real line in three regions, with a neighborhood of the classical turning point in the 
middle. The classical turning point is defined as the point where the potential is equal to the energy 


afc = ^ log fc . 


(80) 


The WKB solution to first order in the classically allowed region ] — oo,afc — e] (with e an appropriately chosen real 
number, see below) is 


ci = + j) . (81) 

and reduces to a plane wave in the allowed region for a <C ctfe. The presence of the barrier fixes the amplitude and 
the phase relation of the incoming and the reflected wave through the matching conditions. 

The solution in the classically forbidden region [ak -f e, -l-oo [ is instead exponentially decreasing as it penetrates the 
potential barrier and reads 


(-,111 _ _ ^ _ 


exp 6"^“ — A:2 4- ^ arctan ^ 




(82) 


Finally, in the intermediate region ]ak — e, -I- e[ any semiclassical method would break down, hence the Schrodinger 
equation must be solved exactly using the linearized potential 


V{a) ~ V{ak) + V'{ak){a - at) = 4e"‘“''(a - ak) = Ak^{a - ak) ■ 


In this intermediate regime, Eq. (22) can be rewritten as the well-known Airy equation 




- tci^ = 0 , 


(83) 


(84) 


where the variable t is defined as 


t = (4A:^)3 (a - ak) ■ 


(85) 


Of the two independent solutions of Eq. (84), only the Airy function Ai(A) satisfies the boundary condition, and hence 
we conclude that 


Ck = c Mt), 


( 86 ) 


where c is a constant that has to be determined by matching the asymptotics of Ai(A) with the WKB approximations 
on both sides of the turning point. Hence we get that, for A <C 0 


while for A ^ 0 


Ai(A) 


cos(ilAr-f) 


Ai(A) 


2v^Ai/4 ■ 


(87) 


( 88 ) 


We thus fix c = Note that the arbitrariness in the choice of e can be solved, e.g. by requiring the point ak — e 

to coincide with the first zero of the Airy function. 

The WKB approximation improves at large values of A:, as one should expect from a method that is semiclassical 
in spirit. Yet it allows one to capture some effects that are genuinely quantum, such as the barrier penetration and 
the tunneling effect. 
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Figure 5: The blue and the yellow curve represent, the WKB approximation of the wave function in the classically 
allowed and in the forbidden region, respectively. They both diverge at the classical turning point. The Airy 
function is represented in red. The value k = 220 was chosen. 


From the approximate solution we have just found, we can construct wave packets as in Ref. m- We thence restrict 
our attention to the classically allowed region and the corresponding approximate solutions, i.e. Cl{a) for a <C Ofc = 
^logfc and compute the integral in Eq. (25). If the Gaussian representing the amplitudes of the monochromatic 
modes is narrow peaked, i.e. its variance 6^ is small enough, we can approximate the amplitude in Eq. (81) with that 
corresponding to the mode with the mean frequency k. Thus, introducing the constant 


1 


{k^ - e4“)i/4 7ri/4y& ’ 


(89) 


we have 


•0(a, (j)) ~ C-, 


f 


dfc 


(k-ky 
2b2 COS 


I - ^/k'^ — e4“ — - arccoth ( 

V2^ 2 V 




,4a 


4 


ik(p 


(90) 


Moreover, 


k arccoth ( , _ I 

Vv^ 


p4“ , 


: k arccosh ~ k arccosh 

g2a g2a 


(91) 


The last approximation in the equation above holds as the derivative of the inverse hyperbolic function turns out to 
be much smaller than unity in the allowed region. Furthermore, the term approximated in Eq. (91) dominates over 


the square root in the argument of the cosine in the integrand in the r.h.s. of Eq. (90), so we can consider the latter 
as a constant. Hence, we can write 


ip{a, (j)) ~ cj 


dfc e ' 2 i, 2 * cos (Afc + (5) 


(92) 
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where we have introduced the notation 


A = - arccosh , 

2 e2“ ’ 



(93) 

(94) 


for convenience. Using Euler’s formula we can express the cosine in the integrand in the r.h.s. of Eq. (92) in terms of 
complex exponentials and evaluate In fact, defining 




i(fe(A+(^)+(5) 


Qk — e 


—i(fe(A-0)+(5) 


(95) 


we have 




Shifting variables and performing the Gaussian integrations we get the result as in Ref. [10] 


'tP{a,(t))~Ck^b(e '■2 (A+0) V_ + e ‘’.(A <t>? . 


(96) 


(97) 


B. ORTHOGONALITY OF THE MACDONALD FUNCTIONS 

It was proved in Refs. [3TVI33] that 



(/cx) 


- y') 
2 v smh{'Kv) 


(98) 


which expresses the orthogonality of the MacDonald functions of imaginary order. Performing a change of variables, 
the above formula can be recast in the form 


da K.h 



Tr'^S{h — h') 
/isinh(7r^) ’ 


(99) 


which is convenient for the applications considered in this work. For large values of k it is equivalent to the normal¬ 
ization used in Section El 
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5{k-k') . 
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